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Abstract 

The abBA collaboration proposes to conduct a comprehensive pro- 
gram of precise measurements of neutron /5-decay coefficients a (the 
correlation between the neutrino momentum and the decay electron 
momentum), b (the electron energy spectral distortion term), A (the 
correlation between the neutron spin and the decay electron momen- 
tum), and B (the correlation between the neutron spin and the decay 
neutrino momentum) at a cold neutron beam facility. We have used a 
GEANT4-based code to simulate the propagation of decay electrons 
and protons in the electromagnetic spectrometer and study the en- 
ergy and timing response of a pair of Silicon detectors. We used these 
results to examine systematic effects and find the uncertainties with 
which the physics parameters a, b, A, and B can be extracted from an 
over-determined experimental data set. 

Key words: Neutron beta decay asymmetry parameters, detector 
Monte Carlo simulation, GEANT4. 



1 Introduction 



The abBA collaboration is proposing to perform a measurement of a "com- 
plete set" of correlations in the neutron /3-decay using the same apparatus, 
and improve the precision of the correlation coefficients a, b, A, and B by 
up to an order of magnitude. 

GEANT4 is a general-purpose software package for simulation of the 
passage of particles through matter that provides a complete set of tools 
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Figure 1: The schematic drawing of the main elements of the abBA spec- 
trometer. 

for all domains of detector simulation [1J. In particular, the GEANT4 
toolkit currently provides particle tracking in non-uniform magnetic and 
electric fields and handles combined electromagnetic fields transparently l2l 
The GEANT4 Low Energy Electromagnetic Physics group validates the 
low energy electromagnetic processes for electrons down to 250 eV [3J. 

In this report we describe a GE ANT4 simulation of the abBA spectrom- 
eter and outline the algorithm for the extraction of the physics decay pa- 
rameters. 



2 abBA Detector Geometry 

In the tentative design of the abBA spectrometer the decay particles (elec- 
trons and protons) are guided by the electric and magnetic fields and in- 
teract only in sensitive detectors thus avoiding the energy losses and scat- 
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terings in apertures, grids, or windows |4"I . 

The simplified geometry of the detector defines two sensitive planar 
Silicon detectors with a 100 x 100 mm 2 area and a 2 mm thickness. The two 
Si detectors are separated by 4 meters. The coordinate system is defined 
with the Si detectors at x 1:2 = ±4 meters (Fig. [l}. 

A passive solenoid magnet is placed around the decay region, with 
its axis of symmetry perpendicular to the incident neutron beam, along 
the x coordinate. The 3 m long magnet with a 0.8 m radius can produce 
a 4 T central magnetic field that decreases to 1 T at the detector positions, 
thus guiding charged particles from the decay region to the Si detectors. 
A tubular electrode held at ~30 kV accelerates the protons so they can be 
detected in a Silicon detector. 



3 Magnetic Field 

The magnetic field along the x axis of the detector solenoid is given by: 

L oc x \ 

^R 2 + (L-xY + y/WTx^ j ' 

where L is the length of the solenoid, R its radius, and x the axial coordi- 
nate. Meanwhile, N denotes the number of turns per unit length, and / is 
the electrical current in the closely wound cylindrical coil. 

Thanks to axial symmetry, the magnetic field off-axis, outside its sources, 
can be represented in terms of the magnetic field strength B(x, 0) along the 
axis: 

«.(*■ o> = - a / x = t w b,2 ""m (f r = b <*> - j b "^ °> + ■ • 

n=0 ^ '' 

and 

' n=l v ' 



s 2nNI 
B(x,p = 0) = 
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where </> is the magnetic scalar potential and p = \/y 2 + z 2 is the axial 
radius coordinate. These fields have been programmed into the GEANT4 
user routine. 



4 Event Generator 

The electrons from the neutron /3-decay are generated from 5x5x5 mm 3 
central volume with the relativistic differential decay rate given by [5|: 

dT = (G F V ud ) 2 F(E e ) \p e \ E v 2 

dE e dtt p JVL Pv ' (2vr) 5 mn [ Ep + E u + E e 0- p v )\ 

where E e and p e {E u and p v ) are the electron (neutrino) energy and momen- 
tum, m n is the neutron mass, Gf is the Fermi constant, V u d is the Cabbibo- 
Kobayashi-Maskawa matrix element, (3 = p e /E e , and F(E e ) is the Fermi 
function that describes the interaction of the electron and the recoil proton. 
The transition matrix element squared \M\ 2 is given by 

|M| 2 = m n m p E e E u (l + |- e«) (l + |- <#>) 

xC (E e )(l + 3~g\) jl + (l + |- C 1 (E e )P -p v + b 

+ (l + £ *< 2) ) [C 2 (E e ) + C 3 (E e )P ■ p v \h ■ 
+[C 4 (E e ) + C 5 (E e )p-p v }n-p„}, 

where m p is the proton mass, a is the fine structure constant, ey is a low 
energy constant, 5 a 's are model-independent radiative corrections, is 
the axial coupling constant, and the correlation coefficients a, A, B are 
incorporated into the recoil corrections Ci(E e ) 0. 
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abBA EXPERIMENT abBA SIMULATION 
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Figure 2: A flow diagram of the analysis. Monte Carlo histograms are 
recalculated in each step of MINUIT minimization. The "experiment" de- 
pends only on measured quantities (subscript E), while the MC column 
depends both on generated variables (subscript T) and simulated detector 
response (subscript E). 
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5 Results and Conclusions 



A GEANT4 simulation of abBA detector energy and timing response was 
performed for 10 6 neutron /3-decays. We used the values of the correlation 
coefficients from Ref. Q (a = -0.1039, b = 0, A = -0.1161, B = 0.9878). 
For each event we recorded the neutron polarization, generated momenta 
of the final state particles and measured energy depositions and timing 
hits in the Silicon detectors. 

The separate GEANT4 run which included systematic effects (energy 
and timing resolutions of Si detectors, detector calibration uncertainties, 
detector response nonlinearities, magnetic field inhomogeneities, neutron 
polarization uncertainty, etc.) was used to simulate the experimental data. 
("MC data"). The flow chart of the physics analysis is summarized in 
Fig. |2| The correlation coefficients and their fitted uncertainties are ex- 
tracted using the standard MINUIT code |5|. 

We present two example results: (i) the extraction of the parameter a 
with unpolarized neutron beam in Fig. |3l and (ii) the asymmetry coeffi- 
cient A for 80 % polarized neutron beam in Fig. |4| At the current stage of 
development, a GEANT4 simulation limited to 10 6 neutron decay events 
and 10 6 MC data events, runs 24 CPU hours on a 1 GHz Linux computer. 
Given the limited event statistics, analysis of MC data results in the coeffi- 
cient b = -0.0025 ± 0.0028 and the coefficient A = -0.1170 ± 0.0010, where 
statistical and systematic uncertainties are combined. The code will be 
made faster by using the adiabatic invariants for charged particle tracking 
in the electromagnetic field, which will markedly improve the uncertain- 
ties of our method. 
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MINUIT FIT: b = -0.0025±0.0028, x 2 =0-2 
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Figure 3: Monte Carlo energy spectrum of neutron /?-decay electron (top) 
and fractional differences A = (exp — the) / (exp + the) between the simu- 
lated and "MC data" spectrum (bottom). 
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Figure 4: Monte Carlo histogram of the j3 asymmetry (top) and fractional 
differences between the simulated and "MC data" spectrum as a function 
of the electron energy (bottom). 
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